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Summary. — 

Quenching of gluonic jets and heavy quark production in Au+Au collisions at RHIC 
can be understood within the pQCD based 3+1 dimensional parton transport model 
BAMPS including pQCD bremsstrahlung 2f>3 processes. Furthermore, the devel- 
opment of conical structures induced by gluonic jets is investigated in a static box 
for the regimes of small and large dissipation. 

PACS 25. 75. -q - Relativistic heavy-ion collisions. 
PACS 25 . 75 . -q - quark-gluon Plasma. 
PACS 24.10.Lx - Monte Carlo simulations. 
PACS 52.35.Tc - Shock waves and discontinuities. 
PACS 25.75.Cj - heavy quark production in HIC. 
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1. Introduction 

The values of the elliptic flow parameter i> 2 measured by the experiments at the 
Relativistic Heavy Ion collider (RHIC) [1, 2, 3] suggest that in the evolving QCD fireball 
a fast local equilibration of quarks and gluon occurs at a very short time scale < 1 
fm/c. This locally thermalized state of matter, the quark gluon plasma (QGP), behaves 
as a nearly perfect fluid, confirmed by viscous hydrodynamics [4, 5] and microscopic 
transport theory [6, 7]. The viscosity to entropy ratio coefficient ry/s has to be rather 
small, possibly close to the conjectured lower bound r]/s = \/Air from a correspondence 
between conformal field theory and string theory in an Anti-dc-Sitter space [8]. To 
achieve a rather small 77/s value in a partonic gas, binary perturbative QCD (pQCD) 
processes require unphysical large cross sections and thus inelastic radiative interactions 
[9] become important. 

The phenomenon of jct-qucnching has been another important discovery at RHIC 
[10]. Hadrons with high transverse momenta are suppressed in Au + Au collisions with 
respect to a scaled p+p reference [11, 12]. This quenching of jets is commonly attributed 
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to energy loss on the partonic level as the hard partons produced in initial interactions 
are bound to traverse the QGP created in the early stages of heavy-ion collisions (HIC). 
In addition, very exciting jet-associated particle correlations have been observed [13], 
which might be the result of a conical emission of propagating shock waves in form of 
Mach Cones induced by highly-energetic partons traversing the expanding medium [14]. 
Furthermore, heavy quarks are interesting for the investigation of the early stage of the 
QGP. Due to their large mass, initially produced heavy quarks, can cover - depending on 
their production point - a long distance through the QGP. Interactions on this way and 
subsequent modifications on heavy quark distributions can reveal valuable information 
about the properties of the medium. 

A large class of phenomena in heavy-ion collisions can be investigated within the 
framework of the kinetic transport theory. Among others, the kinetic transport model 
BAMPS (Boltzmann Approach to Multiparton Scatterings) [9] was developed to describe 
the early QGP of a HIC. Using BAMPS early thermalization of gluons within r < 1 fm/c 
was demonstrated in Au+Au collisions at ^snn = 200 GeV employing Glauber initial 
conditions and the coupling constant a s = 0.3. In addition to the elastic pQCD gg <-> gg 
processes, pQCD-inspired bremsstrahlung gg o ggg was included. This was shown to 
be essential for the achievment of local thermal equilibrium at that short time scale. 
The fast thermalization happens also in a similar way using a Color Glass Condensate 
as inital conditions [15]. 

BAMPS has been applied to simulate elliptic flow and jet quenching at RHIC energies 
[16] for the first time using a consistent and fully pQCD-based microscopic transport 
model to approach both key observables on the partonic level within a common setup. 
The left panel of Fig. 1 shows that the medium simulated in the parton cascade BAMPS 
exhibits a sizable degree of elliptic flow in agreement with experimental findings at RHIC 
as established in [6, 7]. 

The extraction of the shear viscosity over entropy density ratio rj/s has confirmed the 
essential importance of inelastic processes. Within the present description bremsstrahlung 
and back reaction processes lower the shear viscosity to entropy density ratio significantly 
by a factor of 7, compared to the ratio when only elastic collisions are considered [17, 18]. 
For a s = 0.3 one finds ij/s = 0.13, where for as — 0.6 the values matches the lower bound 
of r\js = 1 /4-7T from the AdS/CFT conjecture. 

In this paper we demonstrate the description of different phenomena in relativistic 
HIC using the relativistic pQCD-based on-shcll parton transport model BAMPS. Due 
to the large momentum scales involved the energy loss of partonic jets can be treated 
in terms of perturbative QCD (pQCD) and most theoretical schemes attribute the main 
contribution to partonic energy loss to radiative processes [19]. In addition, the possible 
propagation of Mach Cones in the QGP induced by such highly-energetic partons can be 
studied. Considering the earlier work investigating the effects of dissipation on relativistic 
shock waves [20, 21, 22] we demonstrate the transition of Mach Cones from ideal to the 
viscous one. It is a major challenge to combine jet physics on the one hand and bulk 
evolution on the other hand within a common framework. In the end, the heavy quark 
production at RHIC and the Large Hadron Collider (LHC) during the evolution of the 
QGP is studied. 

2. Jet Quenching in Au+Au collisions at 200 AGeV 

For simulations of Jet Quenching in heavy ion collisions the initial gluon distributions 
are sampled according to a mini-jet model with a lower momentum cut-off po = 1.4 GeV 
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and a if-factor of 2. The test particle method [9] is employed to ensure sufficient statis- 
tics. Quarks are discarded after sampling the initial parton distribution since currently 
a purely gluonic medium is considered. To model the freeze out of the simulated fireball, 
free streaming is applied to regions where the local energy density has dropped below a 
critical energy density e c (e c = 1.0GeV/fm 3 unless otherwise noted). 

The right panel of Fig. 1 shows the gluonic Raa simulated in BAMPS for central, 
b = fm, collisions. It is roughly constant at -R^° ns ~ 0.053 and in reasonable agreement 
with analytic results for the gluonic contribution to the nuclear modification factor Raa 
[19], though the suppression of gluon jets in BAMPS appears to be slightly stronger. We 
expect improved agreement in future studies when employing a carefully averaged (b) 
that will be better suited for comparison to experimental data than the strict b = fm 
case. 




part 



Fig. 1. - Left panel: Elliptic flow V2 as a function of the number of participants for Au+Au at 
200 AGeV for different combinations of the strong coupling a s and the critical energy density 
e c . See [7] for more information. 

Right panel: Gluonic Raa at midrapidity (y e [—0.5, 0.5]) as extracted from simulations for 
central Au+Au collisions at 200 AGeV. For comparison the result from Wicks et al. [19] for the 
gluonic contribution to Raa and experimental results from PHENIX [23] for ir° and STAR [10] 
for charged hadrons are shown. 



We have computed the gluonic Raa for non-central Au + Au collisions at the RHIC 
energy of y/s = 200AGcV with a fixed impact parameter b = 7 fm (Fig. 2), which roughly 
corresponds to (20 — 30)% experimental centrality. A comparison in terms of the magni- 
tude of the jet suppression for b = 7fm is difficult since there are no published analytic 
results available to compare to. Taking the ratio of the b = 7 fm to the b = fm results as a 
rough guess indicates that the decrease in quenching is more pronounced in BAMPS com- 
pared to experimental data. The ratio of the nuclear modification factor between central 
(0 — 10)% and more peripheral (20 — 30)% collisions is Raa\ %_ w % / Raa\ 2 q%-3o% ~ 0-6 
for the experimental data, while for the BAMPS results RAA\t=ofm I ^4^lb=7fm ~ 0-4- 
However, the issue of detailed quantitative comparison needs to be re-investigated once 
light quarks and a fragmentation scheme are included into the simulations. 

To complement the investigations of Raa at a non-zero impact parameter b = 7 GcV, 
we have computed the elliptic flow parameter vi for gluons at the same impact parameter 
and extended the range in transverse momentum up to roughly pt ~ 10 GeV, see left 
panel of Fig. 2. For this a critical energy density e c = 0.6GeV/fm 3 has been used, 
in order to be comparable to previous calculations of the elliptic flow within BAMPS. 
The V2 of high-pT gluons is at first rising with px, but from px ~ 4 to 5 GcV on, it 
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Fig. 2. - Left panel: Elliptic flow V2 for gluons in simulated Au+Au collisions at 200 AGeV with 
b = 7fm. e c = 0.6GeV/fm 3 . 

Right panel: Gluonic Raa as extracted from BAMPS simulations for b = Ofm and b = 7fm, 
e c = 1.0GeV/fm 3 . For comparison experimental results from PHENIX [23] for tt° are shown 
for central (0 - 10)% and off-central (20 - 30)% collisions. 



begins to slightly decreases again. This behavior is in good qualitative agreement with 
recent RHIC data [24] that for charged hadrons shows i>2 to be rising up to t>2 ~ 0.15 at 
Pt ~ 3 GeV followed by a slight decrease. 

3. Transition from ideal to dissipative Mach Cones 

There is an important issue in recent studies of relativistic heavy-ion collisions (HIC) 
whether the small but finite viscosity allows the development of relativistic shocks in form 
of Mach Cones. Within the framework of BAMPS studies were finished to answer the 
question whether shocks can develop with finite viscosity and how this will alter such a 
picture [20] . Within the relativistic Ricmann problem it was shown that one dimensional 
shocks smears out if viscosity is large [21, 22]. However, the expected viscosity in HIC 
seems to be small enough to allow a significant contribution of shocks in form of Mach 
Cones into the picture of HIC. In the following we report a very recent study. 

Mach Cones, which are special phenomena of shock waves, have their origin in ideal 
hydrodynamics. A very weak perturbation in a perfect fluid induces sound waves which 
propagate with the speed of sound c s = y/ dp/de, where p is the pressure and e is the 
energy density. In the case where the perturbation with velocity Vj e t propagates faster 
than the generated sound waves, the sound waves lie on a cone. Considering a gas of 
massless particles, where e = 3p and c s = l/y/3, then the emission angle of the Mach 
Cone is given by a w = arccos(c s /wj e t) = 54, 73°. 

A stronger perturbation induces the propagation of shock waves exceeding the speed 
of sound, therefore the emission angle changes and can be approximated by 
a rj arccos(?; s hockA>jct)- Here 



(1) 



^shock 



(Pi -Po)(eo +Pi) 



(ei - e )(ei +p ) 



is the velocity of the shock front, po (eo) the pressure (energy density) in the shock 
front region and p\ (ei) in the stationary medium itself. Eq.(l) has the following limits: 
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For po >> pi we obtain v s hock ~ 1, whereas for a small perturbation, po ~ Pij we get 

^shock ~ Cs • 

We employ the microscopic transport model BAMPS to investigate Mach Cones with 
different strength of dissipations in the medium using a jet moving in positive z-direction, 
initialized at t = frn/c at the position z = —0.8 fm. The jet is treated as a massless 
particle with zero spatial volume and zero transverse momentum, that is, p z = E- iCt = 200 
GeV and t>j ct = 1. The energy and momentum deposition to the medium is realized via 
collisions with medium particles. In this scenario we neglect the deflection of the jet and 
it can not be stopped by the medium; its energy and momentum is set to its initial value 
after every collision. 

All simulations are realized within a static and uniform medium of massless Boltzmann 
particles and T = 400 MeV. For this study we consider only binary scattering processes 
with an isotropic cross section among the bulk particles. To save computational runtime 
we reduce our problem to two dimensions. Here we choose the xz-plane and apply a 
periodic boundary condition in y-direction. 



T|/S = 1/64JI T)/S = 1/471 T|/s = \l% 




z[fm] 

Fig. 3. - (Color online) Scenario of a massless jet with p z = Ej e t = 200 GeV which can not be 
stopped by the medium - the shape of a Mach Cone shown for different viscosities of the medium, 
rj/s — 1/64-7T (left), rj/s = 1/4tt (middle), rj/s = 1/ir (right). We show the energy density plotted 
together with the velocity profile. Additionally, in the left panel the linear ideal Mach Cone for 
a very weak perturbation is shown by a solid line; its emission angle is a w = 54, 73°. 

In Fig. 3 we demonstrate the transition from ideal Mach Cone to a highly viscous 
one by adjusting the shear viscosity over entropy density ratio in the medium from 
r]/s = 1/647T w 0.005 to 1/tt rj 0.32. The energy deposition of the jet is approximately 
dE/dx = 11 — 14 GeV/fm. We show a snapshot at t = 2.5 fm/c. 

Using an unphysical small viscosity of rj/s = 1/647T we observe a strong collective 
behavior in form of a Mach Cone, as shown in the left panel of Fig. 3. Due to the fact 
that the energy deposition is strong, the shock propagates faster than the speed of sound 
through the medium. For comparison, the ideal Mach Cone caused by a very weak 
perturbation is given by a solid line with its emission angle a w = 54, 73°. Furthermore, 
a strong diffusion wake in direction of the jet, characterized by decreased energy density, 
and a head shock in the front are clearly visible. 

If we increase the viscosity of the medium to larger values, shown in the middle and 
left panel of Fig. 3, the typical Mach Cone structure smears out and vanishes completely. 
Due to stronger dissipation, the collective behavior gets weaker because of less particle 
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Fig. 4. - (Color online) Scenarion of a deflectable jet with finite energy p z — Ej et = 20 GeV - 
Induced Mach Cone structure for different viscosities of the medium, rj/s = 0.01 (left), r)/s = 0.1 
(middle), 77/s = 0.5 (right). We show the energy density plotted together with the velocity 
profile. 



interactions in the medium with a larger 77/s. The results agree qualitatively with earlier 
studies [20, 22] , where a smearing-out of the shock profile is observed with higher viscosity. 

In addition to the scenario of an unstoppable jet we demonstrate in Fig. 4 the scenario 
of a massless jet with finite energy which can also be deflected. Its initial energy is set 
to p z = £jct = 20 GeV, where the starting point is z = —0.3 fm. We show the results for 
different viscosities, 77/s « 0.01 to 0.5 at t = 5.0 fm/c. In analogy to the results above we 
observe a clear Mach Cone structure for small viscosities and a smearing out with larger 
values of 77/s. Only in the ideal case a strongly curved structure in which the building 
up of a strong vortex is visible. The physical meaning of these phenomena and also jets 
with the full pQCD cascade have to be explored in future studies. 

4. Heavy quarks in BAMPS 

Initial heavy quark production during hard parton interactions in nucleon-nucleon 
scatterings and secondary production during the evolution of the QGP are studied. We 
use the event generator PYTHIA [25] to determine the initial heavy quark distributions, 
which agree with the experimental data from PHENIX [26]. Nevertheless, these dis- 
tributions have large uncertainties due to their sensitivity on the parton distribution 
functions in nucleons, the heavy quark masses as well as the renormalization and factor- 
ization scales (see [27] for a detailed analysis). For the initial distribution of the gluonic 
medium we use three different approaches: the mini-jet model, a color glass condensate 
inspired model and also PYTHIA in combination with the Glauber model. 

In the following we give a brief overview of our results on heavy quark production in 
the QGP. More details concerning this section can be found in [27]. 

Secondary heavy quark production in the QGP is studied within a full BAMPS sim- 
ulation of Au+Au collisions at RHIC. According to our calculations the charm quark 
production in the medium lies between 0.3 and 3.4 charm pairs, depending on the model 
of the initial gluon distribution, the charm mass and whether a K = 2 factor for higher 
order corrections of the cross section is employed. However, compared to the initial yield 
these values are of the order of a few percent for the most probable scenarios. Therefore, 
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Fig. 5. - Left panel: elliptic flow V2 of heavy quarks with pseudo-rapidity \rj\ < 0.35 at the end 
of the QGP phase for Au+Au collisions at RHIC with an impact parameter of 6 = 8.2 fm. For 
one curve the cross section of gQ — > gQ is multiplied with a K factor. For comparison, data 
of heavy flavor electrons [28] is also shown. In contrast to the left panel, where the theoretical 
curves for V2 on the quark level are plotted, the right panel depicts the flow of heavy quarks 
and electrons, which stem from the decay of D and B mesons produced by fragmentation. 



one can conclude that charm production at RHIC in the QGP is nearly negligible. 

At LHC, however, the picture looks a bit different: Here the charm production in 
the QGP is a sizeable fraction of the initial yield and is even of the same order for some 
scenarios (with mini-jet initial conditions for gluons with a high energy density). In 
numbers, between 11 and 55 charm pairs are produced in the QGP. 

Bottom production in the QGP, however is very small both at RHIC and LHC and can 
be safely neglected. As a consequence, all bottom quarks at these colliders are produced 
in initial hard parton scatterings. 

The elliptic flow and the nuclear modification factor 

d 2 N AA /dp T dy 

(2) V 2 = ( n - ^ ) , Raa - 




N hia d 2 N pp /dp T dy 

(p x and p y are the momenta in x and y direction in respect to the reaction plane) of 
heavy quarks at mid-rapidity are observablcs which are experimentally measurable and 
reflect the coupling of heavy quarks to the medium. A large elliptic flow comparable to 
that of light partons indicates a strong coupling to the medium. On the other hand a 
small Raa is a sign for a large energy loss of heavy quarks. Experimental results reveal 
that both quantities arc on the same order as the respective values for light particles 
[28, 29, 30]. 

As we have recently shown [31] elastic scatterings of heavy quarks with the gluonic 
medium using a constant coupling a s — 0.3 and the Debye mass for screening the t 
channel cannot reproduce the experimentally measured elliptic flow. In order to explain 
the data one would need a 40 — 50 times larger cross section than the leading order 
one. Of course, this K factor is too large to represent the contribution of higher order 
corrections. However, as we demonstrated in Rcf. [32] and is shown in the following, 
the discrepancy with the data can be lowered - even on the leading order level - by a 
factor of 10 by taking the running of the coupling into account and by improving the 
incorporation of Dcbye screening. The remaining factor of 4 difference could then indeed 
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stem from neglecting higher order effects, which, however, must be checked in a future 
project. 

The following calculations are done analogously to [32, 33, 34]. An effective running 
coupling is obtained from measurements of e + e~ annihilation and non-strange hadronic 
decays of r leptons [33, 35]. 

Since the t channel of the gQ — > gQ cross section is divergent it is screened with 
a mass proportional to the Debye mass, which is calculated by the common definition 
m 2 D = An (1 + Nf/Q) a s (t) T 2 with the running coupling. The proportionality factor k of 
screening mass and Debye mass is mostly set to 1 in the literature without a sophisticated 
reason. However, one can fix this factor to k = 0.2 by comparing the dE/dx of the Born 
cross section with k to the energy loss within the hard thermal loop approach toK« 0.2 
[33, 34]. 

Fig. 6 depicts the Raa of heavy quarks, which shows for K = 4 the same magnitude 
of suppression as the data. 

These improvements lead to an enhanced cross section which also increases the elliptic 
flow. The left panel of Fig. 5 shows vi as a function of the transverse momentum px for 
the leading order cross section without any improvements, with the running coupling, 
with the corrected Debye screening and with both modifications. 

The elliptic flow of the latter reproduces the order of magnitude of the data, if the 
cross section is multiplied with K = 4, which is much smaller than the previous employed 
K = 40 — 50 and lies in a region which could account for higher order corrections. 
However, one has to check if these corrections have indeed a similar effect as a constant 
K factor of 4. Therefore, the calculation of the next-to-leading order cross section is 
planned for the near future and will complement 2 f-> 3 interactions for gluons, which 
are already implemented in BAMPS [9]. The shapes of the theoretical curve and of 
the data points are, however, slightly different. This is an effect of hadronization and 
decay to electrons, which is not shown in the left panel of Fig. 5. If one takes those two 
effects into account, the agreement of the data and the theoretically curve is much better, 
in particular for high pt (see right panel of Fig. 5). We performed the fragmentation 
of charm (bottom) quarks to D (B) mesons via Peterson fragmentation [36] and used 
PYTHIA for the decay to electrons. At low pt the agreement between the theoretical 
curve and data becomes worse, since Peterson fragmentation is not the correct model 
here and coalescence may play a role. 
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